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Knotted and tangled structures frequently appear in physical fields, but so do mechanisms for 
untying them. To understand how this untying works, we simulate the behavior of 1,458 superfluid 
vortex knots of varying complexity and scale in the Gross-Pitaevskii equation. Without exception, 
we find that the knots untie efficiently and completely, and do so within a predictable time range. 
We also observe that the centerline helicity - a measure of knotting and writhing - is partially 
preserved even as the knots untie. Moreover, we find that the topological pathways of untying knots 
have simple descriptions in terms of minimal 2D knot diagrams, and tend to concentrate in states 
along specific maximally chiral pathways. 


Tying a knot has long been a metaphor for creating sta¬ 
bility, and for good reason: untangling even a common 
knotted string requires either scissors or a complicated 
series of moves. This persistence has important conse¬ 
quences for filamentous physical structures like DNA, the 
behavior of which is altered by knots and links mm- An 
analogous effect can be seen in physical fields, e.g., mag¬ 
netic fields in plasmas or vortices in fluid flow; in both 
cases knots never untie in idealized models, giving rise to 
new conserved quantities umi. At the same time, there 
are numerous examples in which forcing real (non-ideal) 
physical systems causes them to become knotted: vor¬ 
tices in classical or superfluid turbulence 00, magnetic 
fields on the surface of the sun [7], and defects in con¬ 
densed matter phases [8]. This presents a conundrum: 
why doesn’t everything get stuck in a tangled web, much 
like headphone cords in a pocket [9]? 

In all of these systems, ‘reconnection events’ allow 
fields to untangle by cutting and splicing together nearby 
lines/structures (Figure [^) [2l[7l[8l [T q 14T2] . As a result, 
the balance of knottedness, and its fundamental role as a 
constraint on the evolution of physical systems, depends 
critically on understanding if and how these mechanisms 
cause knots to untie. 

Previous studies of the evolution of knotted fields have 
been restricted to relatively simple topologies or idealized 
dynamics [aumnnii]. Here, we report on a systematic 
study of the behavior of all prime topologies up to nine 
crossings by simulating isolated quantum vortex knots in 
the Gross-Pitaevskii equation (GPE, eqn.[^. We observe 
that all knots untie, regardless of topological complexity 
or scale, and do so along preferred topological pathways. 
As the vortices untie, the loss in topology is compensated 
by a gain in the ‘coiling’ of the unknotted vortices. Cru¬ 
cially, these results are determined by the geometry of 
the vortex knots, rather than the details of the vortex 
reconnections. 

The quantum counterpart of smoke rings in air, 
vortices in superfluids or superconductors are line-like 
phase defects in the quantum wavefunction, 'ip{x) = 




FIG. 1: Reconnection events and vortex knots, a, A 

schematic of a vortex reconnection event, in this case con¬ 
verting a trefoil knot (K3-1) to a pair of linked rings (L2al). 
b. An ‘ideal’, or minimum rope-length, trefoil knot, c. Using 
the center-line of an ideal knot provides a consistent, uniform 
geometry for any knot or link; nearby strands are exactly 
spaced by the rope diameter, drope, which becomes the char¬ 
acteristic radius, ro, of the loops which compose the knot, d. 
Example ideal configurations of topologies with different min¬ 
imal crossing number, n. The number of topologies excluding 
mirrored pairs is indicated in square brackets, e, A 2D slice 
of the phase field of a superfluid wavefunction with a knotted 
vortex line (light blue), f. Example minimal knot diagrams; 
in each case the topology cannot be represented by a simpler 
planar diagram. The chirality of each crossing is indicated. 
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yj , where p and (j) are the spatially varying den¬ 
sity and phase (Figure [^). This quantum wavefunction 
can be mapped to a classical fluid velocity and density 
via the Madelung transform: u = V0; p = [15]. 

A simple description of the time evolution of this super¬ 
fluid wavefunction is given by the Gross-Pitaevskii equa¬ 
tion m; in a non-dimensional form it is given by: 

( 1 ) 

where in these units the quantized circulation around 
a single vortex line is given by: F = § di • u = 27r. 
The GPE has a characteristic length-scale, known as the 
‘healing length’, which corresponds to the size of the 
density-depleted region around each vortex core (C = 1 
in our non-dimensional units if the background density is 
Po = !)• The GPE is a useful model system for studying 
vortex dynamics: vortex lines are easily identified, recon¬ 
nections occur without divergences in physical quantities, 
and the topological dynamics were recently shown to be 
comparable to real viscous fluids m- Moreover, wave- 
functions in the GPE can be numerically evolved using 
a simple split step method [T0| (see supporting online 
material). 

Due to the difficulty associated with generating initial 
states with vortices of arbitrary shape, previous stud¬ 
ies of knotted superfluid vortices have been restricted to 
particular geometries within one knot family: the torus 
knots m- Here, we directly integrate the flow field of a 
classical fluid vortex to produce phase fields with defects 
(vortices) of any topology or geometry |T7| (Eigure[^). 
Using this technique, we are able to observe the evolu¬ 
tion of all prime knot and links with 9 or less crossings, 
a total of 322 distinct topologies. 

To ensure consistency between different topologies, we 
choose the ‘ideal’ form of each knot, equivalent to the 
shape of the shortest knot tied in a finite thickness rope 
(Eigure[^-d) [T8l[T9]. These canonical shapes are known 
to capture aspects of the knot type as well as approxi¬ 
mating the average properties of random knots [20]. To 
test the sensitivity of the evolution to scale and perturba¬ 
tions, we simulate the dynamics at three different scales: 
ro/C = {15,25,50}, where vq is the characteristic radius 
of loops in the ideal knots (Eigurej^). We also consider 
four randomly perturbed copies of each n < S knot/link 
(4 X 123 configurations with r.m.s. deviation a = 0.25ro 
and ro = 15^). We label topologies using a generalized 
notation following [21], e.g. a “stevedore’s knot” is K6- 
1, with the ‘K’ indicating it is a knot (vs. a link, ‘L’), 
n = 6 is the minimal crossing number (Eigure [^), and 
the remainder indicates an arbitrary ordering. 

Eigure shows the evolution of a 6-crossing knot, 
K6-2, as it unties. The knot can be seen to deform to¬ 
wards a series of vortex reconnections that progressively 
simplify the knot until only unknotted rings (unknots) re¬ 
main. This behavior has previously been observed for a 


handful of simple knots and links, here we find the same 
behavior in every of the 1,458 simulated vortex knots. 
Significantly, this unknotting proceeds regardless of scale 
or distortions, indicating that it is a generic phenomenon. 

To quantify the dynamics of vortex knots as they un¬ 
tie, we compute their length, energy, and centerline he- 
licity as a function of the evolving geometry of the vortex 
lines. The non-dimensional ‘centerline helicity’ - which 
measures the total linking, knotting, and coiling in the 
field - is given by [S] dT] [23] : 

h = + ( 2 ) 

i/j i 

where Lkij is the linking number between vortex lines i 
and j, and Wri is the 3D writhe of line i, which includes 
contributions from knotting as well as helical coils. Sim¬ 
ilarly, the vortex energy, can be computed purely in 
terms of the vortex geometry, up to a logarithmic core¬ 
correction factor [24]. Note that this energy is not the 
total energy of the superfluid, which would also include 
sound waves and is conserved in the GPE. 

Eigure [^-f shows the length and helicity of all 322 
topologies as they untie (see also Eigure SI in the sup¬ 
porting materials). Three general trends can be clearly 
discerned from our results: 1) the timescale for unknot¬ 
ting is determined predominantly by the overall scale of 
the knot, 2) the helicity is not simply dissipated, but 
rather converted from links and knots into helical coils, 
with an efficiency that depends on scale, and 3) the vor¬ 
tex lines stretch by ^20% as they untie, even though 
the vortex energy decreases slightly. Interestingly, all of 
these results appear to be independent of complexity: 
simple knots untie just as quickly as complicated ones, 
and lose the same relative amount of helicity and energy 
(Eigure S2). 

A histogram of the unknotting times, Eigure [^-d, is 
consistent with a log-normal distribution. We find that 
once the time is appropriately rescaled, the mean unknot¬ 
ting times for each simulation group are in the range: 
{tunknot) ~ (3.5 — 4.0) Vq/T {vq is the diameter of the 
‘rope’ in which the ideal knot is tied). 

Despite the fact that the vortex knots untie, they do 
not lose all of their initial centerline helicity. After each 
reconnection event, helices with a range of length scales 
are produced on the reconnected vortices. Without any 
spatial cutoff, this process is expected to exactly conserve 
helicity [mus], however, small helices (compared to the 
healing length) are radiated away as sounds waves. As 
a result, we observe an average helicity loss which has 
an empirical Ah/ho oc (ro/C)~^‘^ trend, consistent with 
previous observations of trefoil knots HU- 

In all cases, the total vortex length increases until 
the knot is completely untied, at which point it approx¬ 
imately stabilizes. By contrast, the vortex energy is 
nearly constant except during reconnections, which dissi¬ 
pate a small amount. As with helicity, the relative energy 
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FIG. 2: Geometric evolution of vortex knots, a, The untying of a randomly distorted 6-crossing knot (K6-2, ro = 50^) to a 
collection of unknotted rings. The rescaled time, t' = txT/ vq, is shown for each step. The top section shows density iso-surfaces 
of the wavefunction (red, |'0|^ = 1/2) and the transparent surfaces (teal or purple) show a constant phase iso-surface. Each 
volume has been centered on the vortex, which would otherwise have a net vertical motion; only 48% of the simulation volume 
is shown, b-d. The relative length, vortex energy, and helicity as a function of time for six different ideal links/knots with 
ro = 25^. Each topological jump is marked with a diamond, and the line color changes to gray after the knot has untied, 
e-f, 2D histograms of relative length and helicity as a function of time for all prime topologies with n < 9. The dashed lines 
indicate average values. The helicity histogram, e,, only includes the 269/322 topologies with ho > 1. 


loss decreases with increasing knot scale, however, it fol¬ 
lows an approximate AEy/Eo oc (ro/0~^‘^ trend and so 
is not proportional to the loss in helicity. 

If one assumes that concentrated vorticity distribu¬ 
tions will always expand, these results have an intuitive 
explanation. Collections of unknotted vortex rings may 
separate without stretching individual vortex lines, but 
this is not possible for a linked or knotted configuration. 
At the same time, if the system is undriven the vortex 
lines must change configuration to conserve energy as 
they increase in length: as previously observed for sim¬ 
ple knots, the formation of closely spaced, anti-parallel 
vortex pairs reduces the energy per unit length m- As 
the stretching continues, these anti-parallel vortex lines 
are driven closer together until they ultimately reconnect; 
this process continues until the knots are completely un¬ 
tied. Interestingly, such a picture also naturally produces 
the anti-parallel reconnection geometry that favors helic¬ 
ity conservation. 

Although the above results demonstrate the over¬ 
whelming tendency for vortex knots to untie, they do not 
elucidate the specific topological pathways which produce 
this untangling. To measure these unknotting sequences, 


we identify the topology, T^, of the vortices after each 
reconnection by computing their HOMFLY-PT polyno¬ 
mials [26l [27] (see supporting materials). This process 
reduces each simulation to an ordered list of visited knot 
types (e.g. Figure]^), allowing us to consider the decay 
process in terms of topological invariants. Due to the 
high degree of symmetry of ideal knots, reconnections 
are often nearly coincident in time, preventing identifi¬ 
cation of the intermediate topology. To avoid this com¬ 
plication, we only consider the decays of the randomly 
distorted knots, which break this symmetry. 

The first question we examine is whether the knot is 
simplifying at each step. We quantify the knot complex¬ 
ity via the crossing number, n, of each knot in a minimal 
2D diagram (Figure El)- Table shows the statistics of 
the jumps in the crossing number through all reconnec¬ 
tions, revealing knots are about an order of magnitude 
more likely to ‘untie’ (An < 0) than ‘retie’ (An > 0) at 
each individual reconnection. On average, more than one 
crossing is removed with each reconnection, underscoring 
the fact that physical reconnections of vortices in 3D are 
not equivalent to removing (or adding) a single crossing 
from a 2D minimal knot diagram. Nonetheless, the min- 
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FIG. 3: Histograms of the rescaled untying time a-d, and the untied vs. initial helicity e-h, for four different groups of 
simulations: a-c, e-g, all 322 ideal knots with n < 9 at a scale of ro = {15, 25, 50}^ and d,h, four randomly distorted versions 
of each n < 8 ideal knot with ro = 15^ and a — 0.25ro (492 simulations total), a-d. The distribution of untying times is 

well described by a log-normal distribution (dashed red line): P{t) oc ^ exp [^— j, where exp /x = {4.0, 3.9, 3.5, 3.7} and 

cr = {0.37, 0.41, 0.44, 0.47} for a-d,, respectively, e-h. The hnal helicity is approximately proportional to the initial helicity (red 
line). The degree to which helicity is preserved depends on overall scale, but is not slightly affected by randomly distorting the 
knots. 


An = -1: 34.4% An =-2: 41.8% An <-3: 12.5% 

An = +1: 3.1% An = +2: 5.8% An > +3: 0.3% 

Removes crossings of same sign (|An| = |Air|): 96.1% 
Jump ends in maximally chiral {\wfinai \ = rifinal)- 82.6% 
Jump leaves maximally chiral: 0.7% 

TABLE I: Probabilities of topological jumps of various types. 

imal diagrams reveal a clear trend towards topological 
simplification. 

If each reconnection does not correspond to ‘removing’ 
a single crossing from a 2D knot diagram, is it still possi¬ 
ble to produce an intuitive description of these events in 
terms of such diagrams? This question can be answered 
by considering the 2D topological writhe, 'ie(T^), which is 
obtained by summing the handedness (±1) of each cross¬ 
ing in a minimal knot diagrams (obtained from [28} 1 ^ ). 
Remarkably, we find that the vast majority (96.1%) of 
reconnection events only remove crossings of the same 
sign from 2D diagrams, i.e. |An| = |Au;|. Such reconnec¬ 
tions have a simple interpretation: they are equivalent 
to the relaxation of a parallel or anti-parallel pair in a 
2D diagram (Figure]^. The removal of multiple cross¬ 
ings thus occurs by a single reconnection happening in an 
anti-parallel pair, followed by the untwisting of a topolog¬ 


ically trivial loop by type-I Reidemeister moves [MIEI]. 
Reconnections followed by more complicated simplifica¬ 
tions are possible, however such events are observed to 
be very rare. 

Figure shows the topological writhe and crossing 
number of every knot with n < 8 (including non-prime 
topologies), connected by lines indicated the frequency 
of the observed unknotting pathways. In addition to il¬ 
lustrating the above results, this diagram reveals the im¬ 
portance of the ‘maximally chiral’ topologies, for which 
\w\ = n. The topological writhe for any particular knot 
or link is bounded by the number of crossings; maximally 
chiral knots and links saturate this bound, which corre¬ 
sponds to every crossing having the same sign. 

Despite the fact that only around a third of all n < 8 
topologies are maximally chiral, 82.6% of jumps end in 
such a state. The dominance of this pathway has a sim¬ 
ple interpretation: if we assume all reconnections satisfy 
I An I > I Are I (corresponding to a slope of |An/An;| > 1 
in Figure [33], once the vortex knot decays into a max¬ 
imally chiral topology it can only leave such a state by 
increasing its crossing number. Indeed, due to the ‘gap’ 
between maximally and non-maximally chiral states, the 
crossing number must increase by An > +2 to leave the 
maximally chiral branch. Moreover, even in the event 
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b relaxation of anti-parallel pair 



1 reconnection + untwisting 
(2 X Reidemeister type I) 


C relaxation of parallel pair 



3 separate reconnections 


FIG. 4: Topological mechanisms for untying vortex knots, a, Knot diagrams of an observed decay pathway for a 6 crossing 
knot; all steps can be described by local untwisting events. The minimal crossing number, n, and topological writhe, ic, is 
labeled for each diagram, b, c. Nearly all reconnection events can be classihed either as relaxation of a twisted pair in an 
anti-parallel or parallel orientation; in either case An = — |Aic|. Reconnections of anti-parallel pairs are equivalent to a crossing 
removal plus one or more type-I Reidemeister moves. 


that the crossing number does increase by this amount, 
we observe that it still typically stays on the maximally 
chiral branch. Thus, statistically, most knots are fun- 
neled into maximally chiral pathway during their unty¬ 
ing, after which they decay only along this pathway. 

Our observation of a preferred maximally chiral path¬ 
way is a generalization of a previously known result for 
site-specific recombination of DNA knots: any p = 2 
torus knot/link (which are all maximally chiral) may only 
convert into another p = 2 torus knot via reconnections 
if the crossing number is decreasing [32]. Our results in¬ 
dicate that this torus knot pathway is one example of a 
more general phenomena. Intuitively, this suggests un¬ 
tangling knots tend to end up in states which are twisted 
in only one chiral direction. 

Taken as a whole, we find that the topological behavior 
of superfluid vortex knots - even complicated tangles - 
can be understood via simple principles. All vortex knots 
untie, and they tend to do so efficiently: monotonically 
decreasing their crossing number until they are a collec¬ 
tion of unknotted vortices. This suggests that non-trivial 
vortex topology in superfluids - or any fluid with similar 
topological dynamics - should only arise from external 
driving. Even in the presence of driving, the observed 
decay pathways indicate that vortices would likely settle 
into a maximally-chiral topology; it would be of great 
interest to probe for such states in superfluid or classical 
turbulence. 

Our results relate the global geometry of vortices and 
their topology, and thus should be independent of the 
small-scale details of superfluid vortex reconnections. 
Fundamentally, analogous reconnection events determine 
the evolution of topology in a wide range of fields, result¬ 
ing in the already noted connections to classical viscous 
fluids and DNA. The geometric and topological nature 


of our results suggests that they might apply more gen¬ 
erally, forming a universal set of mechanisms for under¬ 
standing the role knots play in a variety of physical sys¬ 
tems. 
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